home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DTRCO.z / DTRCO
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDTTTTRRRRCCCCOOOO((((3333FFFF))))                                                            DDDDTTTTRRRRCCCCOOOO((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DTRCO   - DTRCO estimates the condition of a double precision triangular
  10.      matrix.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE DTRCO(T,LDT,N,RCOND,Z,JOB)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      TTTT DOUBLE PRECISION(LDT,N)
  20.         T contains the triangular matrix.  The zero
  21.         elements of the matrix are not referenced, and
  22.         the corresponding elements of the array can be
  23.         used to store other information.
  24.  
  25.      LLLLDDDDTTTT INTEGER
  26.         LDT is the leading dimension of the array T.
  27.  
  28.      NNNN INTEGER
  29.         N is the order of the system.
  30.  
  31.      JJJJOOOOBBBB INTEGER
  32.         = 0         T  is lower triangular.
  33.         = nonzero   T  is upper triangular.  On Return
  34.  
  35.      RRRRCCCCOOOONNNNDDDD DOUBLE PRECISION
  36.         an estimate of the reciprocal condition of  T .
  37.         For the system  T*X = B , relative perturbations
  38.         in  T  and  B  of size  EPSILON  may cause
  39.         relative perturbations in  X  of size  EPSILON/RCOND .
  40.         If  RCOND  is so small that the logical expression
  41.         1.0 + RCOND .EQ. 1.0
  42.         is true, then  T  may be singular to working
  43.         precision.  In particular,  RCOND  is zero  if
  44.         exact singularity is detected or the estimate
  45.         underflows.
  46.  
  47.      ZZZZ DOUBLE PRECISION(N)
  48.         a work vector whose contents are usually unimportant.
  49.         If  T  is close to a singular matrix, then  Z  is
  50.         an approximate null vector in the sense that
  51.         NORM(A*Z) = RCOND*NORM(A)*NORM(Z) .  LINPACK.  This version dated
  52.      08/14/78 .  Cleve Moler, University of New Mexico, Argonne National Lab.
  53.      Subroutines and Functions BLAS DAXPY,DSCAL,DASUM Fortran DABS,DMAX1,DSIGN
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.